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Although faster transport velocities of coUoid-associated actinides, bacteria, and virus than nonreactive 
solutes have been observed in laboratory and field experiments, some questions still need to be answered. To 
accurately determine the relative velocity ( Up J Ut) of ^"Pu and tritium representative of the bulk water, a 
conceptual model of electrostatic interactions coupled with the parabolic water velocity profile in pore 
channels is developed. Based on the expression for Up J Ur derived from this model, we study the effects of 
water flow rates and ionic strengths on the Up J U-j;- Also, the velocity relationship between Pu, tritium and 
Sr^"^ is explored. The results show that Up J Ut increased fairly linearly with decreasing water flow rates; Up J 
Ur declined approximately exponentially with increasing Na"^ concentrations; the charge properties of 
colloid- associated Pu (negative), tritium (neutral) and Sr^"^ (positive) had a close association with their 
transport velocities as L/pu : U-r : Us/'^ - 1.41 : 1 : 0.579. 



Because of long half-life and high radiotoxicity of ''''^Pu, its mobility is a key concern involved in both 
performance assessments of high-level radioactive waste repositories' '^ and remediation of contaminated 
fields such as nuclear test sites''''. Transport velocity of Pu (Up^) is one necessary parameter for both 
evaluating its mobility and then predicting its potential risks to human health. Organic and inorganic colloids 
with the diameter of 1 nm to 1 |im are ubiquitous in underground and vadose zone water^ The Pu(IV) having 
great sorption affinity for the colloids could thus exist in a stable species by strong association with functional 
groups on the coUoidal surfaces'"". Accordingly, L/pu is a direct reflection of mobilization of the colloid- associated 
Pu in the geosphere. 

More recently, coUoid-associated Pu, Am, and Np were observed to have transport velocities 1-5% faster than 
tritium (U-y) in porous media'"'"'. The tritium (tritiated water) is a neutral tracer and its transport may be 
representative of bulk water velocities. In field and column experiments, model colloids such as bacteria""^', 
virus^^"^'', and nanoparticles^^'^'' were also found to travel much faster than anionic ions (e.g. CI" and Br"). This 
velocity enhancement was commonly explained by size exclusion^^ and matrix diffusion^" mechanisms. In brief, 
colloid- associated actinides and/or colloids have relatively low diffusion in pore channels because of colloid-sized 
effects, and hence are excluded from much less mobile water domains near the pore walls where the water 
velocities are lower. The pore channels include both the pore spaces formed adjacent to grain-grain junctions 
and the intra-grain pores larger than the colloid sizes. 

These velocity enhancement phenomena were widely reported. There are, however, stiU some questions that 
remain unanswered. One is how to determine velocity relationship between Pu (or other actinides) and tritium, 
i.e. the velocity enhancement factor Up J U^, since the previously observed relationships were considered as rough 
estimates'''""'. Typically, the ratio of two locations with the highest effluent peaks in breakthrough curves for 
colloidal contaminants and solute tracers was generally specified as the factors''''"'^"'^"*'^". Such locations were 
usually plotted by time or effluent volumes vs. concentrations. The velocity relationship was thus established 
according to the only two single samples collected successively from the effluents'''"™, and then might have large 
deviations from the true factors. The locations of the highest peaks are, in fact, difficult to be targeted sometimes. 
For example, in the case of slightly faster velocities of contaminants than solutes''^'^', the two specified locations 
most often met at the almost same time. For long-pulse experiments, the effluent concentrations in the steady 
limbs of the breakthrough curves exhibited nearly identical peak heights^'''^° '". In such cases, it is hard to judge 
which two single data should be used for indicating their relative velocity relationship. In this study, a simple 
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conceptual model of electrostatic interactions coupled with the para- 
bolic water velocity profile is developed to more accurately determine 
the UpJUy by using all breakthrough data instead of the two single. 

The other question is whether environmentally relevant condi- 
tions, typically including the variations in water flow rate and ionic 
strength, influence the UpJU^- The pore water flow rates functioned 
as a physical mechanism may dominate colloid-associated Pu and 
tritium dispersion'''' '''' that is directly related to their transport path- 
ways in porous media systems''^'^''"'. The lower dispersion of colloidal 
contaminants than that of the solutes^^'^' ''° can lead to an inconsist- 
ent diversity of their pathways. In addition, the variation in this 
physical condition may give rise to the changes in mobile and 
immobile water contents in pore water"" Accordingly, the various 
UpJUi is expected to arise. On the other hand, the ionic strength as a 
chemical mechanism may control the electric double layer repul- 
sion'*'* ''^ between natural colloids as carriers of Pu and pore walls. 
The natural colloids and media, such as soils and granites, usually 
have negatively charged surfaces under natural aquatic conditions. 
This repulsive interaction may hence alter the transport domain of 
colloid-associated Pu in pore channels but not the neutral tritium. It 
is well known that the water velocity has the parabolic distribution 
across the section of the channels. Therefore, the change in UpJUj 
with ionic strengths would possibly occur. Until now, enough atten- 
tion has not been paid to the dependence of UpJU-y on these two 
environmentally relevant conditions. 

Moreover, since the colloidal forces have the potentials to alter the 
transport domain of Pu, the electrostatic attraction between posi- 
tively charged ions (e.g. Sr^*) and media surfaces may regulate the 
transport domain of Sr^* and then the U^^* IU-y- Strontium-90 as one 
fission production from nuclear activities has been large-scale 
released into subsurface environments'"''''^. It is necessary to deter- 
mine accurately its relative transport velocity by our new model 
rather than the two single data method. 

The aims of this work are to develop one conceptual model for 
establishing the velocity relationship of UfJU-Y, to investigate the 
role of water flow rates and ionic strengths in U^JUi, and to obtain 
the relative velocity of l/pu : Ut : Usr'^ ■ The mineral colloids with an 
average diameter {d^ of 460 ± 78 nm were isolated from sandy soils 
collected at the vadose zone of Lop Nor and used for preparing the 
colloid-associated Pu. The sandy soil grains (300 to 700 |.im) are used 
as the porous media packed into the columns. 

Results 

Electrokinetic properties of colloids and media. The (-potentials 
shown in Fig. SI were determined according to the measured 
electrophoretic mobilities. Negative values of the (-potentials are 
an indication of negatively charged surfaces on the soil colloids 
(460 ± 78 nm shown in Fig. S2) and media (elemental compo- 
sition shown in Table SI). The negative surface charge emanates 
from acidic functional groups'"*, dominating the mineral surfaces 
under environmentally relevant pH 8.5 conditions. The charge 
properties of the mineral surfaces became less negative with 
increasing Na^ concentrations. This suggests that the amount of 
variable negative charge decreased while introducing the cationic 
Na* to the suspensions. Large differences in (-potentials of the 
colloid and media surfaces were not observed, due to the fact that 
they were derived from the same in situ sandy soils. The surface areas 
(As) values are reported in Table S2, resulting from Fig. S3. The 
colloid As (69.4 m^/g) is larger than that of the media (15.7 m^/g), 
leading to higher affinity of the colloids for Pu. 

Relative transport velocity of colloid-associated Pu and tritium as 
affected by water flow rates and ionic strengths. The results of Pu 
oxidation state distribution as determined via solvent extraction (Fig. 
S4) show that Pu(IV), both in solution and on the surfaces, was the 
predominant oxidation state (> 95%) in the suspensions, and that 



minor changes of the distribution were not observed in the effluents. 
The high Pu(IV) percentage was presumably due to the reduction of 
initial Pu(V/VI) in the stock solution to Pu(IV) by the redox-sensitive 
Fe and Mn on the colloid surfaces, as also observed by many 
researchers'"'™. The relative transport of Pu and tritium in the 
porous media as affected by the water flow rates (Q) is shown in Fig. 1. 

Normalized effluent concentrations of Pu (C/Cq) were increased 
with the Q, suggesting that hydrodynamic shear acting on the col- 
loid-associated Pu was an important process to promote colloid- 
facilitated transport of Pu. Notably, Pu concentration fronts 
appeared ahead of tritium. Accordingly, Pu(IV) was transported in 
the colloidal species but not the dissolved species through the porous 
media. This was further confirmed by the colloidal faction of Pu (> 
0.995) in the effluents of the peak positions. This fraction is nearly 
identical to the CJCq = 0.996 of the initial Pu suspensions (see Fig. 
S5 at Na+ = 0.002 mol/L). 

Since Pu eluted earlier than the tritium, the consistent collection 
percentages (Cp) of Pu and tritium could not meet at the same loca- 
tions with the effluent volumes (V) but at the two separate and 
Vpu with Vx > ^Pu- The collection percentage is defined as, Cp = 
(cumulative effluent mass at one effluent volume/the total effluent 
mass) X 100%. In the other word, the time for Pu arriving at one 
specified Cp was earlier than that for tritium. The difference in their 
effluent volumes, i.e. Vs = ^ Vpu, designated as the separation 
volumes, may be experimentally determined by using the break- 
through data in Fig. 1. 

Because of high recovery percentages (i?p = 100%) for CI" and 
tritium as inert tracers fed into the columns, the difference in their V 
at two separate locations where their -Rp(s) were overall 50% was used 
to determine their Vg by Krupp et al.^'. In contrast to the inert tracers, 
Pu had significantly lower C/Co (<C 1) and thus extremely small _Rp 
(e.g. the maximal Rp = 0.324% at Q = 4.56 mL/min). Hence we 
determine the between Pu and tritium at the two locations with 
the effluent volumes of Vt(cp=5o%) and Vpu(cp=5o%) but not with the 
Vt(Kp=5o%) and Vpu(j;p=5o%): 

Vs = V't(Cp = 50%) — ^Pu(Cp = 50%) (1) 

The dash lines in Fig. 1 show the Vx(cp=5o%) and Vpu(cp=5o%)- For the 
inert tracers, their total mass (mj) injected into the columns are 
almost equal to their total effluent mass because of Rp ~ 100% 
so that the specified y with the Cp = 50% can yield also _Rp = 50% 
resulted from Rp = (m^ X 50%)/mi = 50%. Therefore, Eq. (1) is 
considered as a common expression suitable for the reactive as well as 
the inert tracers. 

Plutonium exclusion volumes Vex-Pu, i.e. the separation volumes 
per unit mass of the media, are determined as 

^ex-Pu = — (2) 
m 

where m (45.0 g) is the mass of porous media packed into one 
column. 

The water content in the exclusion domain (R — ri domain in Fig. 
S6), i.e. the exclusion water content Oex-Pui is calculated as 



where a and H are the radius of the column (1.0 cm) and the height of 
the media (10.0 cm). 

The velocity relationship between colloid-associated Pu and tri- 
tium can be derived from a conceptual model of electrostatic inter- 
actions coupled with the parabolic water velocity profile (Fig. S6), 
and thus is found to exist as 

t/pu/l/T=l+%^ (4) 

where Up^ and Ux are the average transport velocity of Pu and 
tritium; 0(cm''/cm'') is the measured water content in the column. 
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Effluent volume, l/(mL) Effluent volume, l/(mL) 

Figure 1 | Breakthrough curves for Pu and tritium in the porous media under the conditions of 0.002 mol/LNa* and different water flow rates (Q). The 

dash lines indicate the effluent volumes where the collection percentages (Cp) of Pu (black) and tritium (red) are overall equal to 50%. The 
difference in their effluent volumes (blue line areas), i.e. Vs = Vt(Cp=5o%) ~ Vfu(Cp=5o%)< is the separation volumes {Vs). The initial trace-level 
concentrations of Pu (Q) = 10"' mol/L. The error bars resulted from the standard deviations of Co (three samples) and from the measured relative 
standard deviations (RSD) of C. 



The derivation of Eq. (4) from the conceptual model developed in 
this study is presented in the Supplementary Information. 

The experimental conditions and the determined parameters of 
^S) ^ex-Pu' Oex-Pu' ^rid UpJU-j; are reported in Table 1. 

The separation volumes Vs rnay further express the exclusion 
domain that was totally inaccessible to the colloid-associated Pu, 
ascribed to the repulsive interaction mechanism. In contrast, tritium 
could diffuse into the entire domain in pore channels and therefore 
its transport velocity (Ux) is indicative of bulk pore water velocities. 
As shown in Table 1, the Vex-pu increases with decreasing Q, so that 
the colloid-associated Pu was forced into smaller spaces surrounding 
the center axis where the pore water velocity was highest. This might 



cause a fast transport velocity of Pu relative to the bulk pore water, i.e. 
UpJ L/t > 1 . Furthermore, the relative velocity of UfJ Ut is increased 
approximately linearly from 1.35 to 1.52 with the Q from 4.56 to 
0.472 cmVmin, as shown in Fig. 2a. The higher UfJU-Y can be attrib- 
uted to the continually compressed transport domain closer to the 
center axis. 

The Oex .Pu was related to immobile water (layer), since the exclu- 
sion domain was inaccessible to the colloid-associated Pu. The 
immobile water content for the tritium {Oim.-T> ^s often determined 
by the Two-region Model) should be Oim-T < Oex-pu> due to its mole- 
cule diffusion faster than the Brownian diffusion of colloid- 
associated Pu into the immobile water domain. Gvirtzman^^ 



Table 1 | The experimental conditions and their effects on separation volumes (Vs), exclusion volumes (Vex-Pu)/ exclusion water content 
(^ex-Pu)/ and relative transport velocity of Pu and tritium [Upu/Lh] 



Flow rate, Q (cm^/ min) 


Na+ (mol/L) 


Water content, 0 (cm^/cm^j 


Vs (mL) 


Vex-Pu (cmVg) 


Oex.Pu (cmVcm^) 


UpJUj 


4.56 


0.002= 


0.372 


4.06 


0.0903 


0.129 


1.35 


3.37 


0.002 


0.345 


4.16 


0.0925 


0.133 


1.38 


2.28 


0.002 


0.306 


4.42 


0.0983 


0.141 


1.46 


0.472 


0.002 


0.280 


4.56 


0.101 


0.145 


1.52 


2.28 


0.01 


0.306 


3.86 


0.0857 


0.123 


1.43 


2.28 


0.17 


0.306 


1.92 


0.0428 


0.0613 


1.21 


2.28 


0.96 


0.306 


1.36 


0.0302 


0.0432 


1.14 



"This concentration was produced while preparing Pu suspensions by the introduction of NaOH to adjust the pH to 8.5. 
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Figure 2 | Responses of the relative transport velocity ( Up J Ux) to (a) water flow rate and (b) Na* concentration in pore water. The error bars resulted 
from the measured RSDs of the effluent concentrations ( C) of Pu and 'H, and the experimentally determined water content (6). 



indicated the existence of t'lm-x < ''ex-anion for the anionic ions such 
as and CI", resulting from their exchange model. Evidently, 

''im-T < Oex-Pu was an indication that Pu was forced into the more 
centre streamlines, so that UpJU-y > 1 was observed. This is con- 
sistent with the results determined via our conceptual model. In 
general, the immobile water content tends to increase with the 
decreases in Q and/or Q*^-*->-^^^ as followed by the O^^-Pu in Table 1. In 
comparison, the increasing trend in ()ex-Pu potentially prevailed in 
terms of the slower Brownian diffusion. This thus led to rising Up J 
U-Y with decreasing Q, as shown in Fig. 2a. 

Size exclusion (i.e. hydrodynamic chromatography) thoroughly 
discussed during the past several decades^'-^^'"'^'', seems to also 
account for the velocity enhancement phenomena. This colloid-sized 
effect, involving preferential exclusion of the coUoid-associated Pu 
from smaller, more tortuous pore channels, could cause straighter 
transport pathways for the Pu. Therefore, the relatively fast velocity 
of Pu was a result of the lower extent of its pathway diversity. 
Nevertheless, size exclusion, emphasizing on the transport pathway 
diversity but not the transport domain in pore channels, is suitable 
for examining the relative movement of neutral particles rather than 
the charged species whose transport domain can be altered by surface 
forces. In addition, the colloidal particles and media inevitably 
exhibit the heterogeneous size distributions in natural environments. 
In such cases, it is difficult to establish the velocity relationship by the 
size exclusion mechanism. These aspects have been considered with 
the variable separation volumes in our conceptual model. 

The increase in UfJU-Y with decreasing Q was a result of the 
increased Oex-Pu- Here, Oex-Pu depends on this physical condition. 
Then, what wiU happen to the Oex-Pu if the change in chemical con- 
ditions occurs in pore water? Further transport experiments were 
performed under the constant physical conditions but different 
Na* concentrations in the initial Pu suspensions, as shown in 
Fig. 3. It is noted that the size exclusion mechanism could not explain 
the dependence of UfJUi on the variation in Na^ concentrations. In 
this study, we can determine the UpJU-r by the Eq. (4) and other 
parameters by the equations from (1) to (3). The results are provided 
in Table 1. 

It is observed in Fig. 3 that the normalized effluent concentrations 
of Pu still appeared ahead of tritium. This indicates that the relative 
transport characterized by UfJU-Y > 1 could also take place in spite 
of the constant Q. The Vs is decreased from 4.42 to 1 .36 mL with Na^ 
concentrations (0.002 to 0.96 mol/L) and the corresponding Vcx-Pu 
from 0.0983 to 0.0302 cm^'/g. It is evident that Na* concentrations 
had an effect on the Pu transport domain in pore channels. 



The electrostatic interaction between the two negatively charged 
surfaces of the media and the colloids may be considered as the sum 
of two contributions, including the van der Waals attraction ^vdw 
and the electric double layer repulsion ^edi (The calculations are 
presented in Supplementary Information). As shown in Fig. 4, the 
^vdw is independent of the ionic strength in pore water, whereas the 
repulsion (^edi) markedly decreases with the increase in Na^ con- 
centrations, regarding the same interaction distances. The decreased 
repulsion could make the colloid- associated Pu closer to the pore 
walls, and thus higher Na^ concentrations resulted in lower U^JU-y- 
The l/pu/Ur, which declines approximately exponentially with 
increasing Na* concentrations, is presented in Fig. 2b. On the other 
hand, the relative high repulsion at lower Na* concentrations could 
result in transport of more coUoid-associated Pu through the systems 
(see C/Co in Fig. 3). 

No studies have been carried out to determine the quantitative 
velocity relationship between Pu and tritium under these environ- 
mentally relevant conditions (Q and ionic strength). Incorporation 
of their velocity relationship into the widely used transport mod- 
els^^"^' is necessary to evaluate the environmental risk posed by the 
low- solubility actinides into the subsurface. 

Relative transport velocity of coUoid-associated Pu, tritium, and 
Sr^*. The results by 10 kD ultrafiltration indicate that the Sr(N03)2 
in the Pu suspensions mainly existed as the dissolved species 
(—98.2%), i.e. the positively charged Sr^* ions^"'^'. The relative 
transport of colloid-associated Pu, tritium, and Sr^^ in one porous 
media system are displayed in Fig. 5. As expected, Sr^* appeared 
behind tritium, revealing that the transport velocity of Sr^* was 
slower than that of tritium. Then, their transport velocities may be 
expressed as l/pu > Ut > t^s/*- 

Due to the attractive electrostatic interaction between Sr""* and 
pore walls, Sr""* could be attached onto the walls, and thus be 
excluded from the center domain in pore channel. This domain 
(illustrated by in Fig. S6) reflects the separation volumes 
between Sr^* and tritium, which is opposite to the Vg between col- 
loid-associated Pu and tritium. 

^S = Vsr2+(cp = 50%)~ Vt(Cp = 50%) (5) 

where Vs/*(cp=5o%) and Vt(Cp=5o%) are the effluent volumes of Sr^* 
and tritium. 

The exclusion volumes Vex-s/* and exclusion water content 
Oex-s/* also determined by Eq. (2) and (3), respectively, are reported 



SCIENTIFIC REPORTS | 4:5037 | DOI: 10.1038/srep05037 



4 



0° 0.0012- 

(2 0.0010- 



= 0.0008- 



0.0006- 



(D 



0.0004- 



0.0002- 



0.0000- 



0,00035 

Q 

a. 0.00028 

o 

c 

9 

1 0.00021 



" 0.00014 



■o 0.00007 



E 

5 0.00000 



1.0 0.00035 
0.00030 




M Plutonium 
Tritium 
Na =O.I7mol/L 



1.0 



0.8 



0.6 o 



0.4 



0.2 



10 20 30 40 50 60 70 



0.0 E 



Effluent volume, l/(mL) 



20 30 40 50 
Effluent volume, l/(mL) 

Figure 3 | Breakthrough curves for Pu and tritium in the porous media under the conditions of water flow rate (Q = 2.28 mL/min) and different 
Na* concentrations. The dash lines indicate the efQuent volumes where the collection percentages of Pu (black) and tritium (red) are overall equal 
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Figure 4 | The van der Waals attraction (^vdw) ^nd electric double layer 
repulsion ( ^edi) between the pore wall and the colloid surfaces, as affected 
by Na* concentrations, responded to the variation in interaction 
distances. The possible effects of Na* on the colloid size and hence on the 
CPvdw were not taken into account in the calculations (see the 
Supplementary Information) . 



in Table 2. The relative transport velocity of Us/'^/Ut is established 
as 

U,^./Ur = l-^-^ (6) 

The derivation of Eq. (6) based on the conceptual transport model is 
presented in the Supplementary Information. 

As a result of the electrostatic attraction, Sr^* ions stayed in the 
domain near the pore walls, and the exclusion domain surrounding 
the center axis had the water content of Oex-Si'* = 0.129 cm^'/cm''. 
The transport velocity of Sr^* relative to that of tritium behaved as 
Ust^^/Ut = 0.579. The Sr^* tailing can be attributed to desorption 
(detachment) of the Sr^* previously attached onto the pore walls, and 
then these desorbed Sr^* slowly moved in the domain near the walls. 
Recently, matrix diffusion of the dissolved species, such as Cl~, I~, 
Cs*, and tritium in porous fractures was intensively studied, and the 
consistent results indicated the larger and smaller diffusion for the 
cations and anions, respectively, as compared with the neutral tri- 
tium™ This carries an implication of slower transport velocity of 
the cationic Sr^* than tritium, which is further demonstrated by 
Us/^/U-Y = 0.579 resulted from our conceptual model allowing for 
the electrostatic attraction. 

The electrostatic interactions between contaminants and media 
surfaces may exhibit as repulsion or attraction, depending on the 
charge properties of the contaminants. In this study, the charge 
properties have proven to have a close association with their 
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transport velocities. The observed relationship between the colloid- 
associated Pu, tritium and Sr^* is quantitatively determined as Upu : 
Ut : l/s/+ = 1.41 : 1 : 0.579. 

Discussion 

We developed a new conceptual model of electrostatic interactions 
coupled with the parabolic water velocity profile in pore channels. 
This model is much more suitable for determining the transport 
velocities of the charged contaminants. Some studies were carried 
out by using the expressions for UpJU-Y and L7s/*/Ut derived from 
this model. The results show that, 

(1) The relative velocity of colloid-associated Pu and tritium {UpJ 
Ui) increased fairly linearly with decreasing water flow rates, 
resulting from the continually compressed transport domain of 
Pu close to the center streamlines; 

(2) The UpJU-Y declined approximately exponentially with 
increasing Na* concentrations because of the decreased elec- 
trostatic repulsion between the natural colloids as carriers of Pu 
and the pore walls; 

(3) The charge properties of colloid-associated Pu (negative), tri- 
tium (neutral) and Sr^* (positive) had a close association with 
their relative transport velocities, which was determined as l/pu 
: U-Y : Us/* = 1.41 : 1 : 0.579 under the experimental conditions 
of Q = 2.27 mL/min, water content 0 = 0.306 cm^'/cm', Sr^* = 
0.0065 mol/L, and Na+ = 0.01 mol/L. 



The previously reported transport velocities of the colloid-assoc- 
iated Pu, Am, and Np were ~ 1-5% faster than tritium'*"'", as deter- 
mined by only two single data in breakthrough curves. This rough 
estimate possibly caused lower velocities than the UpJUi (1.14-1.52, 
depending on the experimental conditions in Table 1 ) determined by 
the newly developed model in this study. Our results do fall into the 
ranges observed from the field experiments in the aquifers, such as 
1.25 times (bacteria/Br")^', 1-2 times (bacteria/groundwater)'"', 1.5- 
2.5 times (bacteria/Br^)^''''"''. Because filed experiments of actinides 
remain scarce, these bacteria transport may be analogous to the 
coUoid- associated Pu by following considerations. Both bacteria 
and coUoid- associated Pu had the nano-scale sizes and the negatively 
charged surfaces, and therefore they exhibited similar transport 
mechanisms such as size exclusion and repulsive interaction. 

Methods 

Porous media. Sandy soils collected in northwestern China (Lop Nor) were air-dried, 
sieved with a stainless steel sieve. The mineralogical composition of the soils 
determined by X-ray diffraction (D/MAX-2500) are 50% quartz, 15% anorthose, 1 1 % 
sericite (clay minerals), 6% orthoclase, 5% calcite, 4% dolomite, 4% chlorite (clay 
minerals), 2% amphibole, 1% gypsum, 1% pyrite and 1% undetected minerals. The 
elemental composition determined by X-ray fluorescence spectroscopy (Axios) is 
presented in Table SI. Other characteristics were a cation exchange capacity of 
0.0286 mmol/g and pH of 8.5 ± 0.1 (1:1 water-media). 

Natural colloids. The natural colloids were isolated from the sandy soils. The sieved 
soil sample of 150 g (< 0.3 mm in diameter) was added to a glass beaker (4 L pure 
water), and ultrasonically dispersed for 15 min. The colloid suspensions with Stoke's 



Table 2 | Separation volume Vs, exclusion volume Vex, exclusion content 0ex/ and relative transport velocity between both Pu/tritium and 
tritium/Sr^^ under the conditions of v^ater flow rate Q= 2.27 mL/min, water content 0 = 0.306 cm^/cm^, Sr^+(Co) = 0.0065 mol/L, 
and Na+(Co) = 0.01 mol/L 

Pu Sr^* Relative transport velocity 



Vs(mL) Vex-Pu(mL/g) 


Oex-Pu (cmVcm^) 


Vs (ml) Vex-sr'^ (ml/g) 


Oex-sr^^ (cmVcm^) 


UpJUj 




3.93 0.0873 


0.125 


4.05 0.0899 


0.129 


1.41 


0.579 
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diameter of < 1 were siphoned from the upper suspension into a polypropylene 
vessel and then stored in the refrigerator (4''C), used as the colloidal source materials. 
Because the isolated colloids were derived from the sandy soils, their elemental 
composition was identical to the reported results of the media {see Table Si). 

To determine the mass concentration of the parent suspensions, five 25 mL ali- 
quots of the suspensions were transferred to concave Teflon membranes, respectively, 
and dried at SO^'C by infrared light. The determined concentration was 267 ± 6 mg/L. 
The parent suspensions could be diluted with pure water to obtain the desired colloid 
concentrations. 

Characterization of colloids and media. Both Brunauer Emmett Teller (BET) 
surface areas (As) and pore volumes {Wg) of the porous media {300-700 |im in 
diameter) and the dried colloids were measured with N2 adsorption method at 77 K 
(ASAP 2020, Micromeritics). Adsorption/desorption curves of N2 are shown in Fig. 
S3, and the measured results are given in Table S2. 

The porous media of about 5 g were crushed in an agate mortar. The crushed 
particles larger than 1 0 |j.m in diameter were removed using gravity sedimentation. 
Electrokinetic potentials (Q of the remained particles < 10 |j.m and the colloids, as a 
function of Na"^ (NaCl) concentration from 0.001 to 1.0 mol/L, were determined 
using Nano ZS (Malvern) under the conditions of pH 8.5, and Na"^ — 0.002 mol/L. 
The experimental results are presented in Fig. Si. 

The colloid diameter in the prepared suspensions was measured with the dynamic 
light scatter (Nano ZS, Marlven) under the conditions of Na^ — 0.002 mol/L, pH ~ 
8.5, and colloid concentration — 100 mg/L. The result in Fig. S2 shows that the 
average diameter (dc) was 460 ± 78 nm, and that the colloidal suspensions could 
remain relatively stable. 

Pu suspensions. Small amounts of ^^^Pu stock solution (95 ng/g) in 1.5% HNO3 were 
drop by drop added to the stirred colloid suspensions in Teflon bottles. The atom 
ratios in the '^^Pu stock solution were """Pu/^'^Pu - 0.0346, ^'*'Pu/"^Pu - 0.000355 
and ^^Pu/^^^Pu — 0.0000323. After 10 min, the pH of the suspensions was adjusted to 
8.5 using 0.5 mol/L NaOH. Such prepared suspensions were the Pu suspensions. 
Tritiated water { —800 Bq/g) was added to the Pu suspensions for studying the relative 
transport of Pu and tritium. The Na"^ concentrations were adjusted by adding NaCl to 
the suspensions, which was designed to study the effects of ionic strengths on their 
relative transport. In addition, Sr(N03)2 was added into one aliquot of Pu suspensions 
containing tritium in order to explore the relative transport of Pu, tritium and Sr^"^. 

Prior to transport experiments, additional experiments were carried out in order to 
examine whether Pu introduced in the colloid suspensions could sorb onto the 
natural colloids. One aliquot (1 mL) of the Pu suspensions was transferred to the 
10 kD ultrafiltration tube (Amicon Ultra-4, Millipore) to determine the concentra- 
tions of the dissolved Pu (C^, mol/L) that was the filtrate fraction passing through the 
10 kD membranes. The colloidal Pu {C^, mol/L), i.e. the colloid-associated species, 
which was retained on the membranes, was thus obtained as — Cq — Cj. The total 
concentration of Pu (Co) was measured before —0.5 h of the transport experiments, 
considering that the potential sorption of Pu on the Teflon walls. These experiments 
were performed in triplicate. The dissolved (Cj/Cq) and colloidal fractions (CJCq) in 
the Pu suspensions, as a function of the Na^ concentrations, are presented in Fig. S5. 

The trace-level concentration of ^^^Pu, colloid concentration, and pH in the Pu 
suspensions to be used for the transport experiments were Cq = 1 X 10"^ mol/L, c — 
100 mg/L, and 8.5, respectively. The initial Sr^"^ concentration was 0.0065 mol/L. 

Column transport experiment. The porous media (300-700 |j.m in diameter) were 
loaded into eight polypropylene columns with the dimensions of 2.0 cm in diameter 
and 10 cm in height, respectively. Aqueous solution free of colloids was fed into the 
columns from the top by one peristaltic pump, until the steady water contents 
monitored by weighing the columns were achieved. Because the porous nylon 
membranes were placed on the base of the columns to maintain the capillary pressure, 
the unsaturated conditions used for simulating the vadose zone were present in the 
systems. Subsequently, one pore volumes (Vp) of the colloid suspensions with a 
concentration of c = 100 mg/L was introduced into the columns to displace the 
existing water in the systems. The Pu suspensions (c — 100 mg/L) of one Vp were 
injected drop by drop into the columns, followed by injection of the colloid 
suspensions for four Vp. The effluents were continually collected by a fraction 
collector at the column outlet. The ^^^Pu concentration, tritium radioactivities, and 
Sr^^ concentrations were measured by the ICP-MS (Element, Finnigan MAT), the 
liquid scintillation counter (Wallac 1414), and the ICP-AES (VISTA-MPX), 
respectively. 

The columns from 1 to 4 were employed to find the relationship between the water 
flow rates and the UpJUj. The change in water flow rates from 0.472 to 4.56 mL/min 
was controlled via adjusting the rotating speed of the peristaltic pump. The rela- 
tionship between the ionic strengths from 0.002 to 0.96 mol/L (Na"^) and the UiiJUj-. 
were examined using the column 3 and the columns from 5 to 7. The relative velocity 
of Pu, tritium and Sr^"^ was studied using the column 8. For each transport envir- 
onment, both aqueous solution and colloid suspensions injected into the columns had 
the same pH (8.5), colloid concentration (100 mg/L), and Na^ concentration as the 
Pu suspensions containing tritium (or Sr^"^). 

All the experiments were carried out under atmosphere condition (0.03% CO2, v/v) 
and room temperature (25 ± l°C). Some transport parameters are the followings: 
loading density (pb = 1.56 g/cm^), pore volumes (Vp = 12.3 cm^), and effective 
porosity (e — 0.41). 



Analytical methods. Plutonium oxidation state distribution, both in solution and on 
the colloidal surfaces was determined by a combination of centrifugal ultrafiltration 
(10 kD, Amicon Ultra-4, Millipore) and solvent extraction (TTA and HDEHP, 
Sigma-Aldrich Co.) according to the analytic flow diagram in Fig. S4. 

Plutonium-239 concentrations of all samples were measured by ICP-MS using the 
isotope dilution method, as follows. Plutonium-242 was used as a spike, and the atom 
ratio of ^^^Pu/^*^Pu was 0.001148. After introducing known amounts of ^^^Pu in the 
samples, HNO3 was added to adjust the solution to 8 mol/L HNO3. NaN02 was used 
to adjust Pu valence to Pu(IV) at 90 'C for 15 min, and then the solution cooled at 
room temperature. These solutions were fed into Dowexl X 2 resins (Dow Chemical 
CO.) packed in glass columns with 3 mm X 50 mm length. U and matrix elements 
in the feeding solution were removed via washing the columns using 5 mL of 8 mol/L 
HNO3, 3 mL of 10 mol/L HCl, and 3 mL of 3 mol/L HNO3. The purified Pu was 
eluted from the columns with 1.2 mL of 0.01 mol/L HNO3-O.OI mol/L HF solution, 
and determined with the ICP-MS. 
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